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ABSTRACT 

We present a method for measuring small deviations from axisymmetry of the potential 
of a filled gas disk. The method is based on a higher order harmonic expansion of the 
full velocity field of the disk. This expansion is made by first fitting a tilted-ring 
model to the velocity field of the gas disk and subsequently expanding the velocity 
field along each ring into its harmonic terms. We use epicycle theory to derive equations 
for the harmonic terms in a distorted potential. The phase of each component of the 
distortion can vary with radius. We show that if the potential has a distortion of 
harmonic number to, the velocity field as seen on the sky exhibits an to — 1 and m+ 1 
distortion. As is to be expected, the effects of a global elongation of the halo are similar 
to an m — 2 spiral arm. The main difference is that the phase of the spiral arm can 
vary with radius. Our method allows a measurement of e po t sin2</32, where e po t is the 
elongation of the potential and ip2 is one of the viewing angles. The advantage of this 
method over previous approaches to measure the elongations of disk galaxies is that, 
by using H I data, one can probe the potential at radii beyond the stellar disk, into the 
regime where dark matter is thought to be the dominant dynamical component. The 
method is applied to the spiral galaxies NGC 2403 and NGC 3198 and the harmonic 
terms are measured up to ninth order. 

The residual velocity field of NGC 2403 shows some spiral-like structures. The 
harmonic analysis indicates that the to = 3 term is dominant, with an average value of 
~ 0.02v c . This is consistent with an average ellipticity of the potential of e pot sin 2if2 = 
0.064 ± 0.003, but spiral arms may couple significantly to this result. 

In the harmonic analysis of the kinematics of NGC 3198 the m = 2 and m = 3 
terms are strongest (~ 0.01w c ). The inferred average elongation of the potential is 
0.019 ± 0.003. Since the amplitude of the elongation is coupled to the viewing angles 
and may be influenced by spiral arms, more galaxies should be examined to separate 
these effects from true elongation in a statistical way. 

Key words: Galaxies: individual: NGC 2403, NGC 3198 - galaxies: kinematics and 
dynamics - galaxies: spiral - dark matter 



1 INTRODUCTION 

It is generally accepted that disk galaxies have massive dark 
halos but little is known about their shape. Dark halos were 
generally modeled as being spherical until Binney (1978) 
proposed that the natural shape of dark halos is triaxial 
in order to explain warped disks. Triaxial dark matter ha- 
los occur naturally in cosmological N-body simulations of 
structure formation in the universe. These simulations also 
indicate that there may well be a universal density profile 
for dark matter halos (Navarro, Frenk & White 1996; Cole 
& Lacey 1996), but the exact distribution of halo shapes is 
as yet uncertain (Katz & Gunn 1991; Dubinski & Carlberg 
1991; Warren et al. 1992; Dubinski 1994). The main problem 
is numerical resolution, but the inclusion of baryons is also 



important. If baryons are included, the baryonic disk that 
forms in the centre of the dark matter distribution tends to 
make the halo more oblate (Barnes 1994; Dubinski 1994). 
So the observed distribution of the shapes of dark matter 
halos can be a constraint on scenarios of galaxy formation 
and subsequent galaxy evolution. Measuring the shapes of 
dark halos can be split in two parts: measurement of the 
ratio c/a, i.e., flattening perpendicular to the plane of the 
disk, and measurement of the intermediate to major axis ra- 
tio b/a, i.e., elongation in the plane of the disk. Polar ring 
galaxies provide clues on the axis ratio c/a, by comparing 
the rotation curve of the disk with the rotation of the in- 
clined ring. Unfortunately, the values found for c/a are not 
unique, since the interpretation of the data is model depen- 
dent, but all measurements now seem to be consistent with 
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c/a < 1 (Whitmore, McElroy, Schweizer 1987; Sackett & 
Sparke 1990; Sackett et al. 1994). For our own galaxy, van 
der Marel (1991, and references therein) tried to determine 
the flattening of the galactic halo using counts and kinemat- 
ics of halo stars. He found that the axis ratio of the galactic 
dark halo must satisfy c/a > 0.34. Finally, measurements of 
the flaring of the H I disks of spiral galaxies may give us a 
direct handle on the flattening of the dark halo of normal 
spiral galaxies. Oiling (1995a) developed a method to mea- 
sure c/a using gas flaring, and found a very high flattening 
of c/a = 0.1 — 0.5 for the nearly edge-on galaxy NGC 4244 



(Oiling 1995b). Preliminary results from Sicking (1997), who 
developed a method for measuring the thickness of inclined 
gas disks, indicate 0.2 < c/a < 0.8 for NGC 3198 with a best 
value of 0.5, whereas the measurements for NGC 2403 are 
difficult and give no conclusive value for the flattening of the 
dark halo. From these measurements, it seems that spher- 
ical halos can indeed be excluded, although no good value 
for the typical flattening of a dark halo has been found yet. 

The elongation of the dark halo (b/a) is an even more 
undetermined parameter. A few attempts have been made 
to determine the typical elongation of disk galaxies on either 
photometric or kinematic basis. On the basis of a large sam- 
ple (13,482) of spiral galaxies from the APM survey, Lambas 
et al. (1992) find that a pure oblate model of spiral galaxies 
does not fit the observations and b/a ~ 0.9 fits best, consis- 
tent with the findings of Binney & de Vaucouleurs (1981). 
However, Lambas et al. remark that this effect may also be 
due to the presence of spiral arms and bars in the disks. On 
the basis of a distance-limited sample of 766 spiral galax- 
ies, Fasano et al (1992) also conclude that spiral galaxies 
cannot be oblate, although early type spirals seem to have 
more triaxial disks than their late type spirals. Furthermore, 
their late type spiral group falls into two parts: those with- 
out bars seem to be consistent with oblate shapes, whereas 
those with bars are not. These two effects may suggest that 
a triaxial bulge or bar might be responsible for the observed 
triaxiality, although an intrinsic elongation cannot be ruled 
out by these observations. 

More detailed imaging of a small sample of spiral galax- 
ies can give a better handle on the effects of spiral arms and 
bars. For this purpose, imaging of the old stellar population 
in the K'(2.2/x) band was done by Rix & Zaritsky (1995). 
They observed a sample of 18 spiral galaxies and found a 
variety of asymmetries in this sample. They estimated the 
potential in the plane of the disk to have an average elliptic- 
ity of 0.045±g;g|, but remarked that this number may only 
be an upper limit to the true ellipticity, since spiral structure 
still couples significantly to their estimated ellipticity. 

A second way to examine elongation of the potential 
of disk galaxies, is to look at the kinematics of the disks. 
For our own galaxy, Kuijken & Tremaine (1994) find on the 
basis of solar neighbourhood kinematics plus kinematics of 
a more global kind (Hi tangent point velocities, Cepheids, 
etc) that our sun lies near the minor axis of the slightly 
elongated (b/a ~ 0.9) Milky Way disk. Franx & de Zeeuw 
(1992) used the scatter in the Tully- Fisher relation to place 
limits on the elongation of spiral disks, since elongation will 
increase the scatter. Based on the observed scatter alone, 
they find that b/a > 0.9, but they argue that it is unlikely 
that the scatter is caused solely by elongation and therefore 
the value is on average likely to be b/a > 0.94. 



The only direct measurement of the elongation of the 
potential of a disk galaxy so far comes from the SO galaxy 
1C 2006 (Franx, van Gorkom & de Zeeuw 1994, hereafter 
FvGdZ), which contains a large Hi ring in the plane of the 
disk. Using information on both the kinematics and the ge- 
ometry of the ring, FvGdZ were able to measure the elonga- 
tion of the potential at the place of the ring. The measured 
elongation is consistent with zero. FvGdZ used epicyclic the- 
ory to determine the geometry and velocity field of a ring 
in a mildly perturbed potential. If the ring is in an elon- 
gated potential, the velocity variation along the ring will 
not be pure circular rotation, but higher harmonic terms 
will be superposed on it. Making a harmonic expansion of 
the velocity field of the ring and interpreting the measured 
harmonics within the framework of epicycle theory, given 
the ring geometry, allows the elongation of the potential at 
the position of the ring to be found. In this paper, we extend 
the method used in FvGdZ for a single ring to the case of a 
filled gas disk in order to study individual halo elongations 
for a larger sample of spiral galaxies. 

With a single ring, one has the advantage that the ring 
can be considered as a single orbit. Therefore, one knows for 
each point in the ring to which orbit it belongs. We will lose 
this advantage in the case of a filled gas disk. Since galactic 
disks are not obviously made out of separate rings, we will 
decompose the Hi disks of spiral galaxies into individual 
rings (the so-called tilted-ring method, Begeman 1987) and 
perform a harmonic expansion on each of these rings (Binney 
1978; Teuben 1991; FvGdZ). The harmonic expansion can 
be done to arbitrary order and it produces amplitudes and 
phases for each separate order. We will investigate if we 
can qualitatively distinguish between the effects of spiral 
arms and of elongation on the velocity field. We will apply 
this method to two spiral galaxies with extended Hi disks, 
NGC 2403 and NGC 3198. From the measurements of these 
galaxies, the strengths and limitations of our method will 
become clear. We will see that in both cases contributions 
of spiral arms to the velocity field are important, but that 
nevertheless limits to the halo elongation can be set. 

The outline of this paper is as follows: In Section 2 
we lay down the theoretical framework for connecting the 
measured kinematic terms to perturbations of the potential. 
We describe the resulting procedure to fit observed velocity 
fields in Section 3. In Section 4 we present some examples of 
the theory, to see what effects may play a role in interpreting 
the measured harmonics. Next, we examine two spiral galax- 
ies: NGC 2403 in Section 5 and NGC 3198 in Section 6. We 
conclude with a discussion on the results from these galaxies. 
In the Appendix, we give a detailed derivation of the equa- 
tions presented in Section 2. Furthermore, we derive what 
the effect of misfitting ring parameters (inclination, position 
angle and centre) is on the measured harmonics. 



2 DESCRIPTION OF A PERTURBED 
VELOCITY FIELD 

2.1 General perturbations 

If the cold H I gas in a disk galaxy orbits in the potential of 
an axisymmetric dark halo, the resulting density field and 
velocity field are axisymmetric as well. This means that the 
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Figure 1. Geometry of the projected orbit. The line (8 ^ s ,i) is 
the line-of sight. The inclination i is measured from the z-axis, 
the axis perpendicular to the plane of the orbit. The angle £> b s 
is measured from the line 8 = and the azimuthal angle tp is 
measured in the plane of the orbit and is zero at the line of nodes, 
which is indicated by the dashed line. 



velocity field will only show pure circular rotation. Observed 
velocity fields of disk galaxies generally can indeed be fitted 
quite accurately by circular motion. Many different effects 
can cause the remaining small residuals. Here, we will anal- 
yse the case of a small distortion in the potential V , which 
can be written as a sum of harmonic components, each of 
which may have a different pattern speed. This is not the 
most general description possible, as it does not include non- 
stationary phenomena such as spiral arms. But it is a good 
description for the motion of gas in a triaxial halo, since the 
potential distortion, i.e., the elongation of the halo, can be 
considered stationary. Furthermore, we will assume that the 
gas moves in a symmetry plane on stable closed orbits. This 
assumption cannot be correct in detail, since this description 
of the gas motion does not allow for shocks in the gas and 
therefore cannot give an accurate description of spiral arms. 
On the other hand, this is a correct assumption if we are 
i nterested in the motion of gas in a n early circular potential 
(Tohline, Simonson & Caldwell, 1982), like the potential of a 



mildly elongated halo. Although these assumptions restrict 
the validity of our description, it will help us to gain some 
understanding of the qualitative characteristics of velocity 
fields. 

The assumptions that the potential contains a small 
stationary distortion and that the gas moves on stable closed 
orbits allow us to use epicycle theory to analyse the velocity 
fields that originate from such a perturbed potential. We will 
solve the equations of motion in this perturbed potential for 
the possible closed orbits and calculate the velocity field of 
these orbits as it would be seen by an external observer. In 
this way we can directly relate the observed velocity field to 
the potential of disk galaxies. 

Let R and 9 be polar coordinates in the rest frame of 
the galaxy under consideration and consider a rotating, non- 
axisymmetric potential in that frame. The potential can be 
written as: 

V{R,9,t)=V {R)+^V m {R)cos{m[e-n p>m t-<t> m (R)]}.(l) 

m 

Here, Vo{R) represents the unperturbed potential and 



Vm(R) cos{m[0— Q Pt mt— <t> m (R)]} is the m-th harmonic com- 
ponent of the perturbation, assumed to be small compared 
to Vo(-R). The phase of the perturbation as a function of ra- 
dius is denoted as <p m (R). Each of the harmonic components 
is allowed to have its own pattern speed, 0. p , m - If the pertur- 
bation has non-zero pattern speed, or <f> m (R) is not constant 
as a function of radius, then the choice of the line 9 — is 
arbitrary. We focus on the velocity field as generated by a 
potential perturbed by a single harmonic component. The 
general case, equation ([j]), can be obtained by adding the 
results for individual harmonic components, since coupling 
terms between different harmonic components are of second 
order. For a single harmonic component, the potential in a 
frame that corotates with the potential perturbation can be 
written as 



V{R,4>) = V (R) + V m (R) cos{rr# - </> m (R)]}, 



(2) 



where <p = 9 — Q. p>m t. The procedure we follow to find the 
closed orbits in this potential is analogous to the procedure 
followed by Binney & Tremaine (1987, eq. 3-107 and fur- 
ther), where they derive orbits in a potential with arbitrary 
m, but without the radially dependent phase we introduce 
here. We solve the equations of motion for a potential per- 
turbed by a single harmonic component in Appendix Al and 



find the set of possible closed loop orbits for any value of m. 
Let us denote the solutions of the equations of motion in the 
unperturbed potential Vo(R) that represent circular orbits 
by Rq and <po. If we then perturb the potential with a single 
harmonic component as indicated above, we find the set of 
closed loop orbits, equation (A.3): 



R{<t>,t) 



Ro 



2 aim COS T] m — 02mi 



1 J?m) 



<p(R,t) = 0a (*) + a^r (aim + 03m) sin 



m - «4m COS Tj m , 



(3) 



where T] m = rj m (R,t) = m[(j>(t) - </> m (R)] « m[</> (t) - 
<j> m {R)] = m{[Q (R) - tlp, m ]t - <j> m (R)} and n (R) is the 
usual circu lar f requency. The amplitudes a xm are defined in 
equations (A4) of the Appendix. 

The resulting velocity field in the non-rotating frame 
for a general rotation curve is 



I'B 



mv c (R)[l - cj(R)] (faim sin?7 m - a 2m cos??™) , 

v c (R) {l + |[(1 - Lu)a am + (a - tu)a lm ] cos7? m (4) 

+ [m(l — Lu)a,4 m — (1 — a)a 2m ] sinr? m } , 

where u>(R) = Q Ptm /Qo(R) and v c (R) = RQ,o(R), the cir- 
cular velocity. We have defined a = d\n[v c (R)]/ din R. This 
velocity field is observed from a direction (9 b s ,i), where i 
is the inclination of the disk and ^obs is defined as the angle 
between the line 9 — and the observer, see Figure 1. 

Introduce ip — 9 — 9 ^ s + 7t/2 = (p — o b s + tt/2, where 
4>ob& is the angle in the rotating frame that corresponds to 
r? b s . The angle ip is measured along the orbit and is zero on 
the line of nodes. It is sometimes called the azimuthal angle. 
It is defined fully by the orientation of the galaxy on the sky 
and it is independent of the internal coordinate system. It 
is straightforward to derive that the velocity field of a disk 
in pure circular rotation is given by tji os = v c (R) sin i cos ip. 

Let us express the line-of-sight velocity as v\ OE = 
J"] c„ cos nip + s n sin nip. In Appendix Al we show that the 
line-of-sight velocity field has the following form: 



V\os = 



[ci cosip + s m _i sin(m— 1)^> + c m _i cos(m — l)tp 
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where 

Cl 
Sm — 1 



+s m+ i sm(ro+l)*0 + Cm+i cos(m+l)^] • 



(5) 



Cm-l = 



S m +1 = 



««(—■ |{[m— (m+l)o; m + a]a im 
+(l-w m )a 3m }sinmfi„ 1 +i {m(l-uj m )a irn 
+[m(l— Lo m ) — 1 + a]a2m} cosmi^ m ), 
v»(|{[m— (m-rT)u; m + a] aim 
+ (l-aj m )a.3 m }cosm<£ m +| {m(l — u) m )a,4 m 
+[m(l—u m ) - 1 + a] a2m } sin mtp m ), 
v*{\{[m—(m—l)u m - a] am, (6) 
-(1 — a; m )a3 m }sinm^ m + i {m(l — cj m )a 4 m 

— [m(l— w m ) + 1 — a]a2m} cosmy> m ), 
c m+ i = u„(— j{[m— (m — l)w m — a]ai m 

— (1— W m )a3m} COS m<(9 m + i {m(l— w)tt4m 

— [m(l— w) + 1 - a] a 2m } sin mtp m ), 

and we have written u» = u c sini and ip m = o b s — 7r/2 — 
4>m(R)- Now we have expressed all observable parameters in 
terms of internal parameters for a potential perturbed by a 
single harmonic term. 

From equation (||) we can conclude that if the potential 
has a perturbation of harmonic number m, the line-of-sight 
velocity field contains onm-1 and m + 1 harmonic term. 
Qualitatively, this conclusion was also inferred by Canzian 
(1993). Furthermore, we can see from the definitions of the 
a xm in the Appendix, that the amplitudes of the m — 1 and 
m + 1 terms behave in the following way: if Oo > fib, the 
m — 1 term is larger in amplitude than the m + 1 term. If 
flo < fib, the situation is reversed. 

In general, we do not a priori know the inclination i (or 
q = cosi), position angle F and centre (x,y) of the galaxy. 
Here, the position angle F is defined as the angle taken in 
anti-clockwise direction between the north direction on the 
sky and the direction of the major axis of the receding half 
of the galaxy. Usually, it is assumed that the velocity field is 
circular, and the values of the inclination and position angle 
are derived from the best fit. If the true velocity field con- 
tains non-circular terms this procedure will produce wrong 
values of i and F and of the position of the centre. If the 
velocity field is expanded on the basis of these incorrect 
viewing angles and centre, the harmonic terms will change. 
In Appendix A2 we derive the effect of errors in these pa- 
rameters: 5q = 5 cosi, ST and centre 8x,Sy. Under these 
incorrect parameters the line-of-sight velocity becomes: 

J Sx Sq 

vi 03 = v*(R) i(l+a) — -(1-a)— costp 



ST 
4 



(3q+~)-a(--q) 



" Sx 
sin ?/S— ( 1 — a) — — cos 2ij) 



(T) 



-(1-a) ^pj- sin 2%j) + (1 - a) ^ cos 3^ 
2qR Aq 

— ( 1 — a)^-(— — 0) sin 3-0 
' 4 q 

+ S Ci(R) cos iip + Si(R) sin iip. 

i>0 



where ip the angle in the plane of the ring that is zero on 



the apparent major axis of the ring and the Si,d etc. are 
defined as in equation (Q). 

The harmonics arising from an m = 2 term will mix 
with the terms caused by a misfit of the viewing angles. 
Since the tilted-ring algorithm has no way to disentangle 
the effects of a misfit of the viewing angles from that of a 
physical m = 2 term, there will be first order differences be- 
tween the best fitting and true ring parameters for m = 2, 
see Appendix (A3. 2). The same is true for an m = 1 term 
and the centre, Appendix (A3.1). Also, an m — 4 term will 
cause C3,S3 terms, as shown in equation (|^). Due to these 
terms, the position angle and inclination will change as well 
according to equation (H), creating new ci,si,C3,S3 terms. 
This affects the measurement of the elongation of the po- 
tential and is an extra source of noise on our measurement. 
In the same way, an m — 3 term can cause similar effects 
as an m = 1 term. For m > 4, the harmonics in equation 
(^) do not change to first order due to small errors in the 
inferred viewing angles and centre. 



2.2 Global elongation 

We now turn to the case of a globally elongated potential, 
as caused by for instance a triaxial halo. In this case we 
have m = 2, u ~ and (j>2(R) ~ const (see also Appendix 
A3.2). We will assume a flat rotation curve (a = 0), since 



rotation curves in the outer parts of galaxies where the dark 
matter halo dominates are basically flat. Furthermore, as 
rotation curve shapes approach solid body rotation (a = 1), 
the errors in the tilted ring fits become large (see Appendix 
3.3) and a reliable determination of the elongation is no 
longer possible. In the flat rotation curve case, the relation 
between the measured ai2 and the ellipticity of the potential 
e pot is ai2 = e P ot and 032 = 2e poi (FvGdZ). 

The measured line-of-sight velocity along each ring in 
the velocity field will have the following form, as derived 
by FvGdZ (note that their equations for c\ and si contain 
typesetting errors): 

(8) 



^ios = ci cos {[> + li sin ip + S3 sin 3^, 



with 



(•1 



si 



+ e po t cos 2y? 2 ], 



2(1 + 2q 2 + 5q 4 



-e po t sin 2(^2, 



0, 

(l-g 2 )(l + 3g 2 ) . 

v* — ; z e DO t sin 2ir>2 , 

2(1 + 2g 2 + 5g 4 ) P r ' 



(9) 



From the 2i, Si and S3 terms we can derive e pot sin 2<p2- This 
is the quantity that we will be looking for in our measure- 
ments of spiral galaxies. 

Unfortunately, we cannot measure ipz separately (the 
cos 2<£>2 term cannot be measured, since we do not know the 
exact value of u*), so all we can measure is the combina- 
tion of viewing angle and ellipticity. Global ellipticity will 
thus result in si and S3 terms that are constant with radius 
(and thus e po t sin 2<^2 will be constant as a function of ra- 
dius as well) , whereas spiral arms will cause si and S3 terms 
that change sign as a function of radius and e po t sin 2(^2 
will wiggle, as we shall see in Section ^. Radial variation in 
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e po t sin 2ifi2 may help us to distinguish between the effects 
of spiral arms and systematic elongation. 

Equation (^) shows that the C3 term should always be 
zero. If the inclination is not equal to the best fitting in- 
clination, the C3 terms will not be equal to zero. Therefore, 
monitoring C3 in the fitting procedure will be a useful way to 
check the correctness of the fitted inclination. Furthermore, 
spiral arms will also influence the measurements of the in- 
clination. As shown in the Appendix, if <p2(R) is constant, 
the fitted inclination for each ring is identical to the true 
inclination for that ring. But, if tp2(R) varies as a function 
of radius, the fitted inclination will deviate from the true in- 
clination. This effect will not be large (typically one or two 
degrees), but it will be measurable. 



3 FITTING VELOCITY FIELDS 
3.1 The method 

In order to determine the variation of different harmonic 
components as a function of radius, we want to fit the ve- 
locity field of a galaxy along individual, independent rings. 
Within the GIPSY package (Groningen Image Processing 
SYstem, van der Hulst et al, 1992) there is a routine that 
fits a set of so-called tilted-rings to the velocity field of a 
galaxy. The code fits a circular model to the velocity field 
by adjusting the ring parameters (kinematic centre (xo,yo), 
inclination i, position angle P and systemic velocity v sys ), so 
we have a list of ring parameters as a function of radius. All 
the ring parameters are simultaneously fitted with a gen- 
eral least-squares fitting routine. In the tilted-ring fit, we 
should take a uniform weighting function since points near 
the minor axis contain information about ellipticity. It is 
best to keep the centre fixed in the fit because the galaxy 
may contain physical C2 , S2 harmonics in its velocity field. 
In that case, a free centre would drift in such a way as to 
make these terms disappear. If a galaxy would be truly lop- 
sided in its kinematics (i.e., the kinematic centre drifts as a 
function of radius), this would show itself in the C2 and S2 
terms, as can be seen from equation (A26). 

After the tilted ring decomposition we make a harmonic 
expansion of the radial velocity along each ring (using the 
ring's individual F, i, xq, yo): 



Potential perturbati 



■!'lo 



u ByB + £„ cos nip + s„ sin nip, 



(10) 



where k is the order of the fit. The coefficients c n ,s n are 
determined by making a least-squares fit with a basis 

{l, cos ip, sin ip, ■ ■ ■ , cos kip, sin kip\ , (11) 

on the data points in each ring. 

We assume here that we can identify each pixel in the 
velocity map with a unique position in the galaxy. Strong 
warps in the outer H I -layers for instance tend to invalidate 
this assumption, since the value of the velocity field in a 
certain pixel may be composed of two or more rings (the 
line-of-sight crosses 2 or more rings) and the velocity field 
will not be accurate there. Mild warping where the viewing 
angles change by only a few degrees will in general not be 
a problem, as long as individual rings do not overlap each 
other. 




Figure 2. Potential perturbations and decomposed velocity 
fields. From top to bottom: the m = 1,2,3 and 4 terms, from 
left to right : a. The potential perturbation, b. the residual ve- 
locity field (total velocity field minus circular velocity), showing 
mainly the m — 1 component. No tilted-ring fit had to be made 
to find the circular velocity, since the circular velocity is a known 
input paramenter in these models, c. The residual field minus the 
m — 1 component, revealing the m + 1 component. To find this 
component a harmonic fit has been made. 

Normally, when creating residual fields, one would also subtract 
the kinematic m = term. In that case, the kinematic m = 
term in the top row, middle column, would be zero and one 
would only see the kinematic m = 2 term. If the centre would 
be a free parameter in the tilted-ring fit, any kinematic m = 2 
component in the residual field would be zero as well. Finally, the 
alternating behaviour of the m = 3 kinematic term (second row, 
third column) is due to the fact that in the tilted-ring fit the 
inclination is a free parameter. Therefore, the C3 term is fitted to 
zero and only the S3 term is visible. 

All plots are on the same scale. The potential ranges from 
—80 km 2 s -2 to 80 km 2 s~ 2 , all velocity components from 
—2 km s _1 to 2 km s" 1 . White denotes positive, black negative 
values. 



3.2 Error estimates 

The errors given by the tilted-ring decomposition code and 
by the harmonic fitting procedure are formal errors of the 
fits and do not take correlation between points due to beam 
smearing into account. Sicking (1997) has derived that the 
formal errors should be multiplied by a factor 



4-K<J X <Ty 

SxSy 



(12) 



where Sx, Sy are the dimensions of a pixel in arcsec and 
a x ,a y are the dispersions of the beam in arcsec. All mea- 
surement errors presented in this paper are corrected for this 
correlation effect. 



4 SOME EXAMPLE VELOCITY FIELDS 

In order to clarify the equations presented in Section 2 and 
in Appendix A, we present some example velocity fields in 
Figure 2. The different terms range from m = 1 to m = 4. 
For the phase of the fields, ip m (R), we chose a logarithmic 
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Table 1. Parameters used for model fields 

quantity symbol value 



inclination 


i 


60° 


obs. angle 


</>obs 


45° 


pattern speed 




km s — 1 kpc — 1 


amplitude 


V m 


80 km 2 s 2 


circular velocity 


Vc 


150 km s" 1 


maximum radius 


Rmax 


18 kpc 


spiral scale length 


Rs 


10 kpc 


pitch angle 


p 


10° 



spiral, ip m (R) = ln(7?/_R a )/tanp. The amplitude of the per- 
turbation has been taken constant throughout the field. In 
Table [l], we summarise the parameters used to create the 
model fields. The fields are then created using equations ^ 
and Q). Subsequently, these line-of-sight velocity fields were 
fitted as described in Section 3, revealing the two harmonic 
components that were hidden in them. 

The example velocity fields are shown in Figure 2. We 
display from left to right 

(i) The potential containing the m term perturbation. 

(ii) The residual velocity field caused by this potential 
perturbation (i.e., the m — 1 plus m + 1 terms together). 

(iii) The m + 1 component of the velocity field. 

We see that indeed an m term in the density causes 
m — 1 and m + 1 terms in the kinematics of the gas. Since 
the pattern speed fi p was taken zero for all models, the 
m — 1 term will be dominant, see Section 2. The m — 2 
potential perturbation demonstrates the effect of fitting a 
circular model to a non-circular velocity field. The fitting 
procedure adapts the inclination until C3 = so that only 
the S3 term remains. In our example, the S3 term varies 
as a function of radius. If the field would have contained a 
global intrinsic ellipticity, m = 2, <P2(R) ~ const, the si and 
S3 terms would have been constant as a function of radius. 
Therefore, the measured S3 would have been constant and 
non-zero as well, but no radially alternating pattern for S3 
would have been visible. 

Since in the fitting procedure the centre was kept fixed, 
we also see a two-armed spiral in the m — 1 case. If we 
would have taken the centre as a free parameter in the fit, 
it would have drifted in such a way as to make the C2 and 
S2 terms disappear. 

We shall encounter the effects discussed above again 
while interpreting the data of NGC 2403 and NGC 3198. 



5 NGC 2403 

5.1 Data description 

NGC 2403 is a nearby (3.25 Mpc Begeman 1987) Sc(s)III 



galaxy ( Bandage & Tammann 1981 ). It ha s been observed 



4 x 12 hours with the WSRT ( picking 1997j ). The data have 
been smoothed to a circular beam of 13 arcseconds. With 
this beam and a pixel size of 5 x 5.5 arcseconds, the /3-factor 
of equation ([l2|) is 3.73. 

NGC 2403 is well suited for our analysis: 



(i) It has a large extent on the sky (Holmberg dimensions 
29' .0 x 15'. 0) 

(ii) It has a favourable inclination (i » 62°) 

(iii) The data is of high signal/noise and has a small beam 

(iv) The galaxy shows no obvious warp 

Figure 3 shows the surface density map of the Hi gas 
in the upper left corner. All data with a S/N < 3 is clipped 
away, in order to have a reliable velocity field. The ellipse 
is the outermost ring (R — 858 arcsec) fitted to the veloc- 
ity and surface density fields. Outside this radius, the data 
has suffered too much from clipping to make a sensible har- 
monic decomposition possible. The centre of all the rings 
is indicated by a red cross and was kept fixed during the 
tilted-ring fit. 

Next to this map we present the deprojected H I surface 
density map. We took an inclination of 62° to deproject the 
galaxy, the average value found by the tilted-ring fits. NGC 
2403 seems to be a bit elongated in the direction of the 
(apparent) minor axis. Some traces of spiral arms can be 
seen. 

At the bottom left the velocity field is shown. It is rel- 
atively regular apart from some small wiggles in the iso- 
velocity contours, probably caused by spiral arms. Next to 
the velocity map is the residual velocity field, created by 
subtracting the fitted circular velocity from the true veloc- 
ity field. The residuals are not large (< 10 km s _1 in absolute 
magnitude) and appear to show some systematic, three-fold, 
structure. The beam shape is indicated in the lower right 
corner. 



5.2 Results of the harmonic expansion 

The first step in our procedure is to expand the velocity field 
of NGC 2403 into individual concentric tilted-rings. The pa- 
rameters i, F and v sya are measured locally for each ring, in 
order to be able to use the exact formalism of Appendix 
A3. We took the width of each ring to be 26 arcseconds, 
i.e., 2 beams. Therefore, the individual points in the plots of 
Figure 4 can be considered independent. Figure 4 presents 
the inclination, position angle and systemic velocity of each 
ring as a function of radius. The fluctuations in the inclina- 
tion i and position angle V are most probably caused by an 
m — 2 spiral arm in the density. In the case of (p2 (R) 7^ 0, 
the measured inclination will in general not be equal to the 
true inclination, as mentioned before. This will cause the 
measured inclination to wiggle around its average value. We 
expect this average value to be very close to the true in- 
clination. The fluctuations in the systemic velocity, Co, are 
of comparable magnitude as the C2,S2 terms, as it should 
be according to Appendix 3.1 if these variations are caused 
by an m — 1 term in the surface density. It should resem- 
ble — £2, but if C2 is also influenced by an m — 3 term in 
the density, this need not exactly be the case. The weighted 
average values of the ring parameters are: 

(i) position angle: 124? 2 ± 0?04 

(ii) inclination: 61?8±0?18 

(iii) systemic velocity: 133.5 ± 0.09 km s _1 

In Figure 4 we also present the result of the harmonic 
expansion of the tilted-rings. The arrow just beyond the last 
measured point denotes the Holmberg radius. This means 
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Figure 3. WSRT observations of NGC 2403. Top left: Map of the Hi surface density. Some spiral-like structure is visible. The ellipse is 
the outermost ring fitted to the data in the tilted— ring decomposition. Top right: The deprojected Hi image of NGC 2403 (deprojected 
with constant inclination). Bottom left: The velocity field. Systemic velocity is at 133.5 km s _1 . The contours range from 20 km s _1 to 
260 km s _1 and are spaced 30 km s _1 apart. Bottom right: The residual velocity field. Here we have subtracted the result of a tilted-ring 
fit from the original velocity field. Some structure can be seen. Note the small amplitude of the residuals (less than 10%) with respect to 
the original field. The beam size is indicated in the box in the lower right corner. 
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Figure 4. Harmonic expansion and ring parameters of NGC 2403 as a function of radius. Average values of the ring parameters are 
indicated by a dashed line. The points are separated by 26 arcsec (2 beams) and are therefore independent. The Holmberg radius is 
indicated by an arrow at the end of the box (at 870 arcsec). The amplitude of the S3 term is the largest with an amplitude « 3 km s _1 . 
This is about 2% of the rotation velocity ci. At the bottom the average amplitude (j4; = 5^ .(cj(j) 2 + Si(j) 2 ) 1 / 2 /j , A\ = y\ s\(j)/j) 
of the harmonic terms is plotted. 



that our Hi observations are completely within the optical 
disk of NGC 2403. 

After an initial rise, the rotation curve of NGC 2403 is 
relatively flat within errors. The £3 term has been fitted to 
zero by the tilted-ring fit, indicating that the inclination has 
been fitted correctly to the rings (but note again that in the 
case of an m — 2 spiral arm, this inclination is not identical 
to the true inclination of the ring). There is some £4 at radii 



smaller than 400 arcsec that may be caused by an m = 3 
component in the potential. This m = 3 component may also 
influence the £2 and S2 harmonics. These S2 and £2 terms 
are relatively large and are probably caused by an m = 1 
component in the potential. The S3 term is the strongest 
harmonic found in NGC 2403. This term is negative at all 
radii. We can calculate e pot sin 2y>2 from si and S3 at each 
radius as (cf. eq. [9]) 
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(i) It has large dimensions on the sky (Holmberg dimen- 
sions 11'. 9 x 4'.9) 

(ii) The inclination is i » 71° 

(iii) The data is of high S/N and has a small beam 

(iv) NGC 3198 has no warp of significance 

In Figure 6 we present the data for this galaxy. The data 
is of somewhat better quality than the data for NGC 2403, 
since the disk of NGC 3198 is completely filled, although 
the beam is larger. The surface density field shows a ring- 
like structure of a few arcminutes, with 3 spiral-like struc- 
tures extending from it. The ellipses and the circle again 
indicate the outermost data used for the harmonic analysis 
(450 arcsec in this case). The cross indicates the position 
of the centre of the rings. The velocity field of NGC 3198 
is quite regular and the amplitudes of the residual field are 
not very large, typically less than 10 km s _1 in both the 
positive and negative direction. Even though these residu- 
als are small, they do show some systematic structure. It is 
worth noting that the residuals are strongest near the centre 
of NGC 3198. 



Figure 5. The elongation measurement as calculated from the 
s\,s-j, and inclination (eq. |l3[ ). It varies as a function of radius, 
but is negative throughout the galaxy. A possible explanation for 
this behaviour is that the average value of e po t sin 2ip2 is caused 
by some global elongation and the wiggling is caused by other 
effects like spiral arms. The average value is —0.064 ± 0.003 



e P otsin2v? 2 = («3 — h) 



(l + 2g 2 + 5q 4 
£i(l-9 4 ) 



(13) 



Figure 5 shows e pot sin 2^2 as a function of radius. It 
shows significant variations, but is negative at all radii. As 
we saw in the numerical examples of the previous paragraph, 
spiral arms would cause a radially alternating pattern of pos- 
itive and negative values for the S3 kinematic term, whereas 
a global elongation would cause the S3 term to be constant 
and non-zero throughout the disk. A possible explanation 
for the behaviour of the measured li,S3 is therefore that the 
average value of e pot sin 2<fi2 is caused by a global elongation 
of the disk and that the wiggles superposed on it are due 
to spiral arms. The weighted average value of e po t sin 2<f2 is 
—0.064 ±0.003. Since e pot sin2iy92 has only physical meaning 
in the case of a more or less flat rotation curve, it should 
only be trusted outside about 200 arcsec. Harmonics higher 
than m = 5 have small amplitudes. 



6 NGC 3198 

6.1 Data description 

The second relatively nearby (9.4 Mpc, Be geman 1987) spi- 
ral galaxy w e will examine is the Sc(rs)I-II (Bandage & Tam- 



marm 1981) galaxy NGC 3198 . This galaxy was observed 
4 x 12 hours with the WSRT ( picking 199"^ . The observa- 
tions were subsequently reprojected to a circular beam with 
FWHM of 18 arcsec. The (3 factor, equation @, for NGC 
3198 is 4.6. 

NGC 3198 also meets our criteria: 



6.2 Results of the harmonic expansion 

We fitted a set of concentric rings to the velocity field of 
NGC 3198. We present the results for the ring parameters 
in Figure 7. The systemic velocity, Co, varies with an ampli- 
tude of a few km s" 1 , about the same amplitude variation 
as the strongest terms in the harmonic decomposition. The 
inclination of the rings, i, is constant within errors. The po- 
sition angle F varies a few degrees. The weighted values of 
the ring parameters found are: 

(i) position angle: 216?0 ± 0?05 

(ii) inclination: 70?5±0?20 

(iii) systemic velocity: 659.6 ± 0.15 km s _1 

Figure 7 also presents the result of the subsequent harmonic 
fit. After an initial rise, the rotation curve is flat within 
errors up to the last measured point at 450 arcsec. This is 
not the outermost limit of the observed H I distribution (see 
Figure 6) but since the disk shows big gaps at radii larger 
than 450 arcsec, the harmonic expansion would become too 
unreliable. The arrows denote the Holmberg radius, so our 
measurements go out to about 1.3 Holmberg radii. 

As predicted by the theory, the £3 term has been fitted 
to zero by the tilted-ring fit. Together with the fact that the 
inclination is constant as a function of radius, this tells us 
that NGC 3198 does not contain strong m = 2 spiral arms, 
although a global elongation of the disk is still possible. The 
£2, S2 terms are quite strong in the inner part of NGC 3198. 
The centre has been fixed at such a value as to minimize the 
£2, S2 terms. As a consequence the I2, £2 terms are relatively 
small in the outer part of NGC 3198. If the centre is left 
as a free parameter in the tilted-ring fit, it drifts by a few 
hundred parsec in order to make the 62, S2 terms zero. In 
other words, NGC 3198 is somewhat lop-sided. 

In the surface density map we saw that there appear to 
be three spiral-like structures in the Hi . This would cause 
some non-zero m = 2 and m — 4 terms in the kinematic 
expansion of the velocity field. Indeed, we see that there is 
some power in the £4, §4 terms that may be caused by an 
m — 3 component in the Hi field. 
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Figure 6. WSRT observations of NGC 3198. Top left: Map of the Hi surface density. Some spiral-like structure is visible outside an Hi 
ring-like structure. The ellipse is the outermost ring that has been fitted to the data. Top right: The deprojected Hi image of NGC 3198 
(deprojected with constant inclination). The structure and Hi ring is better visible here. Bottom left: The velocity field. The systemic 
velocity is at 659.6 km s — 1 . The contours range from 520 km s — 1 to 790 km s — 1 and are spaced 30 km s — 1 apart. Bottom right: The 
residual velocity field. Here we have subtracted the result of a tilted— ring fit from the original velocity field. There is some four-fold 
structure visible. The beam size is indicated in the box in the lower right corner. 
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Figure 7. Harmonic decomposition and ring parameters of NGC 3198 as a function of radius. Average values are indicated as a dashed 
line. The points are seperated by 1.7 beams and are therefore independent. The Holmberg radius is indicated by an arrow. The total 
amplitude, bottom graph, decreases as a function of harmonic number. The systemic velocity cq varies a few km s -1 over the galaxy, 
comparable in amplitude with the largest harmonics. The position angle T varies a few degrees as a function of radius. The inclination i 
is constant within measurement errors. 



Finally, there is a strong S3 term, together with a less can be due to the variation of e pot with radius, or <P2(R) is 
strong si term. From these two terms we deduce e pot sin2y 2 not completely zero, or a combination of both, 
as a function of radius, as plotted in Figure 8. The low 
amplitude of e po t sin 2y>2 is striking: its weighted average is 
—0.019 ± 0.003. This low value may be caused by e pot itself 

being small, or the viewing angle (p2 being such as to make 7 SUMMARY AND DISCUSSION 
sin 2<fi2 much smaller than unity. The variation in e pot sin 2952 

We have extended the FvGdZ formalism for measuring 
elongation of potentials to the case of a slightly non- 
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Figure 8. Elongation e po tsin2^2, as measured from the si, S3 
terms and inclination. Its radial variation is not very large within 
errors, indicating that global ellipticity might be the cause of the 
non-zero measurement, and its average amplitude is low (— 0.019± 
0.003). The arrow indicates the Holmberg radius. 



axisymmetric filled gas disk, which may contain spiral-like 
perturbations. All the details of this derivation can be found 
in the Appendix. This analysis assumes a stationary per- 
turbation in the frame rotating with the perturbation, and 
closed stable orbits. According to this analysis, a perturba- 
tion in the potential of harmonic number m causes an m — 1 
and m + 1 perturbation in the radial velocity field. The m— 1 
component will be dominant in amplitude if fio > &b, oth- 
erwise the rn + 1 component will be dominant. Applicability 
to non-linear phenomena like spiral arms is limited. But a 
small global elongation of the overall potential, as is the case 
with a triaxial dark matter halo, can be analysed with our 
method. In the case of a triaxial potential, we have m = 2, 
small pattern speed and no radial dependence of the phase 
of the m = 2 perturbation. If we assume a flat rotation 
curve, we can deduce e pot sin2(p2- Here, e po t is equal to the 
ellipticity of the potential and tp2 is the angle between the 
minor axis of the elliptical orbit plus a phase term and the 
viewing angle. The term e pot cos 2<fi2 cannot be determined. 
In order to measure this e pot sin 2<p2 and other perturbations 
of velocity fields, we need to expand the gas disk into indi- 
vidual rings. This was done by fitting a tilted-ring model to 
the velocity fields. A tilted-ring model tries to fit a circu- 
lar model to the velocity field. Given the rings found by the 
tilted-ring fit, we made a harmonic expansion of the velocity 
field along each ring. In that way we obtain the harmonic 
terms c„, s n as a function of R, where we took n < 9 for our 
fits. Given Si(-R) and S3 (J?), we can calculate e po tsin2y>2 
as a function of radius. Other perturbations of the velocity 
field, like spiral arms, do not affect the elongation measure- 
ment, unless they contain some m = 2 or m = 4 component 
in the potential. In that case, the global elongation and the 
other perturbing effect cannot be uniquely disentangled. We 



propose that the absolute mean value of e po t sin 2</?2 is an up- 
per limit of the true elongation and the wiggling is caused 
by spiral arms. Since the method relies on the fact that the 
velocity in each point in the velocity field is uniquely de- 
termined, the method cannot be applied in those parts of 
galaxies where strong warping is evident. 

Two spiral galaxies with regular H I velocity fields were 
analysed with this method: NGC 2403 and NGC 3198. Both 
galaxies show less power in the higher order (m > 5) terms 
than in the lower order terms. 

NGC 2403 clearly shows spiral-like behaviour. Its S3 
term is large and negative at all radii and the amplitude 
changes as a function of radius. Two effects can contribute 
to this measurement: a global elongation of the disk causing 
a constant non-zero S3 and spiral arms, as are visible in the 
surface density map, that cause the wiggles. The radially 
averaged |(e pot sin 2^2}! is 0.064 ± 0.003. 

NGC 3198 shows signs of spiral arms in the m — 4 kine- 
matic mode (in the Hi surface density map a three-armed 
spiral- like structure can be seen) . The relatively large §2 , 62 
terms in the central parts may be caused by lop-sidedness 
of the galaxy and the si, S3 terms are most probably caused 
by global elongation. The deduced value for |(e po t sin 2^2) | 
is 0.019 ± 0.003. 

We conclude that both galaxies may show signs of small 
global elongation (< 0.1), although it is hard to separate the 
effects of spiral arms, global elongation and viewing angles. 
On the other hand, the sum of the effects of spiral arms 
and global elongation is small, which means that the two 
galaxies are close to axisymmetry. But in order to quantify 
this result, more galaxies need to be analyzed in this way. 
Work along these lines is in progress (Schoenmakers et al., 
in prep.). 



ACKNOWLEDGMENTS 

It is a pleasure to thank Tjeerd van Albada for stimulating 
discussions, Konrad Kuijken for a careful reading of an ear- 
lier version of the manuscript, and the referee, Cedric Lacey, 
for his insightful comments that significantly improved the 
paper. F.J. Sicking kindly made available the data for both 
NGC 2403 and NGC 3198. 



REFERENCES 

Barnes, J., 1994, in The Formation and Evolution of Galaxies, 
eds C. Munoz-Tunon & F. Sanchez (Cambridge Univ. Press), 
p. 401 

Begeman K., PhD-thesis, University of Groningen, 1987 
Binney J., 1978, MNRAS, 183, 779 

Binncy J., Tremaine S., Galactic Dynamics , 1987, Princeton Uni- 
versity Press 

Binney J., de Vaucouleurs G., 1981, MNRAS, 194, 679 

Canzian B., 1993, Ap.J., 414, 487 

Cole S., Lacey C, MNRAS, 281, 716 

Dubinski J., 1994, Ap.J., 431, 617 

Dubinski J., Carlberg R.G, 1991, Ap.J., 378, 496 

Fasano G., Amico P., Bcrtola F., Vio R., Zeilinger W.W., 1992, 

MNRAS, 262, 109 
Franx M.,dc Zccuw P.T., 1992, Ap.J., 392, L47 
Franx M., van Gorkom J.M., dc Zecuw P.T., 1994, Ap.J., 436, 

642 (FvGdZ) 



© 1997 RAS, MNRAS 000, 000-000 



Measuring non-axisymmetry in spiral galaxies 13 



Gradshteyn I.S., Ryzhik I.M., 1965, Tabic of Integrals, Series and 
Products, New York: Academic 

Holmberg E., 1958, Medd. Lund Obs. Ser. II, 136 

van dcr Hulst J.M., Terlouw J. P., Begeman K., Zwitser W., Roclf- 
sema PR., 1992, "The Groningen Image Processing System, 
GIPSY" , in: Astronomical Data Analysis Software and Sys- 
tems I, (eds. D. M. Worall, C Biemesderfer and J. Barnes), 
ASP Conf. series no. 25, p. 131. 

Katz N., Gunn J.E., 1991, Ap.J., 377, 365 

Kuijken K., Tremainc S., 1994, Ap.J., 178, 421 

Lambas D.G., Maddox S.J., Lovcday J., 1992, MNRAS, 258, 404 

van der Marcl R., 1991, MNRAS, 248, 515 

Navarro J.F., Frcnk C.S., White S.D.M., 1996, Ap.J., 462, 563 

Oiling R.P., 1995a, AJ, 110, 591 

Oiling R.P., 1995b, BAAS, 187, 48.05 

Rix H-W., Zaritsky D., 1995, Ap.J., 447, 82 

Sackett P.D., Rix H-W.,Jarvis B.J., Freeman K.C., 1994, Ap.J., 
436, 629 

Sackett P.D., Sparkc L.S., 1990, Ap.J., 361, 408 

Sandage A., Tammann G.A., 1981, Revised Shapley-Ames Cata- 
logue of Bright Galaxies (Carnegie Inst, of Washington Pub. 
No. 635) 

Sicking F.J., PhD-thesis, University of Groningen, 1997 
Teuben P., 1991, in Warped Disks and Inclined Rings around 

Galaxies, (eds. S. Casertano, P.Sackett, F. Briggs), Cambridge 

Univ. Press. 

Tohline J.E., Simonson G.F., Caldwell N., 1982, Ap.J., 252, 92 
Warren M.S., Quinn P.J., Salmon J.K, Zurck W.H., 1992, Ap.J., 
399, 405 

Whitmore B.C., McElroy D., Schweizer F., 1987, Ap.J., 314, 439 



© 1997 RAS, MNRAS 000, 000-000 



14 R.H.M. Schoenmakers et al. 



APPENDIX A: HARMONIC ANALYSIS OF THE VELOCITY FIELD OF A FILLED GAS DISK 
Al Deprojection of the line-of-sight velocity 

Let R and 9 be polar coordinates in the rest frame of the galaxy and consider a rotating, non-axisymmetric potential in that 
frame. It can be written as: 



V(R, 0, t) = V (R) + V m{R) cos{m[6 - n p , m t 



.(*)]}. 



(Al) 



The potential of a single perturbation in a frame that rotates with the potential perturbation: 

V^-Rorbit, 0orbit) = Vo(-Rorbit) + V m (R olhit ) COs{m[0 O rbit - 0m (iJorbit )] } , (A2) 

where Vo(Rorbit) is the unperturbed potential and tm(-Rorbit) cos{77l[0orbit — 0m(Rorbit)]} is the perturbation. The subscript 
"orbit" denotes that these coordinates describe the point on the closed loop orbit with guiding centre on (R!o,0o)- 

The derivation given in this Appendix is also valid for a summation of perturbations and the final result is a linear addition 
of the results for the individual perturbations. There is some coupling between different terms that occurs when calculating 
Vcf>, equation (A5), but this coupling occurs in second order equations only and can be ignored in first order calculations. 

Following the treatment of Binney & Tremaine (1987, p. 146), we find that the solution for closed loop orbits is: 



Rorbit = Ro + SR ■■ 

0orbit = 00 + 50 = 



Ro (l — ^rp 1 cos?7 — asm sin 77) 



f>o + 



("lm+"3m) 



Sin T) — a,4 m COS Tj, 



(A3) 

Ro and 0o denote those solutions to the unperturbed problem that represent a circular orbit (R = Ro, = 4>o = (Oo — Q P ,m)t). 
Oo is the circular frequency: Qo = y/Vo/R and O p , m is the pattern speed of the particular perturbation. SR and 5cj> denote 
the perturbations on the circular orbit. The extra terms with respect to the analysis of Binney & Tremaine arise from the 
additional phase term m (Rorbit) in the potential. Furthermore: 



77 = 77(771, -Rorbit) 



dim — 



0.2m = 



ra[(fl 
2 



t 



R(){Qo — fip.m) 



.(-Rorbit)] = m[(pO - 0m (.Rorbit)], 



mVm,4>' m (Rorbit) 



(A4) 



fl 4m — 



V m [(2 + m 2 )fi§ + 2(1 - m 2 )ft fi P ,m + m 2 Q. 2 p 



2 

Ao 

rbit ) 



-Ro(fio 



O "l 2 



,m ) x r ! 
\ T7Z. 77 T y m. 



(Oo — Op, m ) 



Ao(f2o — Op,m) 

where differentiation with respect to R is denoted as ', Ao = Ro [kq — m 2 (Qo — Op.m) 2 ] , and Kq = 3f2 2 + ^0", the usual 
epicycle frequency. These equations break down near resonances. 

The velocities in a frame at rest with respect to the observer are: 

v R = mv c (l - u m ) (|ai m sin 77 - a 2m cos 77) , 

tty = v c {l + \ [(1 — w m )a3m — aj m ai m ] cos 77 - [m(l — w m )a4m — a 2m ] sin 77} , (A5) 
and v c (Ro) = Ro^o, the circular velocity. 



where u) m = ■ 

So far we only looked at what happens to a single orbit. A complete velocity field consists of many orbits and the circular 
velocity may be a function of radius. Suppose we are observing the point (R, 0) in the velocity field. What is the orbit that 
goes through this point, in oth er words, what is the (Ro,0o) that corresponds to (R, 0)? 

Let us write equation (A3) as 



Rorbit(t, Ro, 0o) = Ro [1 + r(R , O )] , 0orbit(t| Ro, 0o) = 00 + p{Ro, 0o)- 

Then to first order the coordinates for the guiding centre of the orbit going through (R, 
R 



Ro = 
Then 

v c (Ro 



l + r(ii,0) 
) = v c {R)(l-r 



= <p-p(R,c 
dln« c (R), 



(A6) 
(A7) 



u(Ro) = w(R)(l-r 



dlnR 
dlnw(R), 



dlnR ' 

Inserting these equations in equation 



(A8) 



Ao 



and retaining only first order terms, we find 



v R = mv c (R) [1 —ct) m (R)] (I ai m sin 77 — a^m cos 77) , 
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«</> = v c(R) {l + — w)a3m + (a — w m )oi m ] cosr] — [m(l — w m )a4m — (1 — a)a2m] sin 77} , (A9) 

where a = d ■ The definition of a is such that the relation between the epicyclic frequency ko and the circular frequency 

Oo can be written as kq — 2(1 + a)Q. 2 . 

We now project the velocity field on the sky. The line-of-sight velocity field is 

Vios(-R) = [vr cos(8 - 8 obs ) - sin(# — 6 ohs )] sini = [vr cos((p - </> bs) - v<j, sin(<p — </>obs)] sin i, (A10) 

where <p obB is the angle between the line <f> — and the observer. The angle i is the inclination of the plane of the disk with 
respect to the observer. Introduce the variable ip = 8 — 8 b B + ?=(/> — (pobs + f • This coordinate is zero along the line of 
nodes. Then: 

Vios(-R) = [vr sin?/) + Vtf, cos ip] sini. (All) 

Now replace <p — » ip + (j> b s — 5 in the expressions for vr and ity and expand the line-of-sight velocity in multiple angles of ip. 
Furthermore, define u„ = v c (R) sini and assume to first order <f>o w 0. After the expansion we find : 

Vios = ci cost/" + s m _i sin(m — 1)^ + c m _i cos(m-l)ip + s m+1 sin(m + l)V> + c m +i cos(m + l)i/>, (A12) 
with 

ci = «» 

s m _i = v*(— i{[m — (m+l)u m +a]ai m + (1— w m )a3 m } sinmip m +§ {m(l-aj m )a4m + [m(l— w m ) — 1 + a] a 2 m} cos mip m ), 
c m _i = v*(|{[m— (m+l)w m + Q]a im + (1— u) m )az m } cosrrvp m +\ {m(l~uj m )a,4, m + [m(l— w m ) — 1 + a]a2m}sinm<^ m ), 
s m+ i = v*(|{[m— (m— l)w m — a]ai m — (l—w TO )a3m} sin m^m+ 2 {m(l— w m )a 4m - [m(l— w m )+l - Q]o 2m }cosm^ m ), (A13) 
c m+ i = v*(— i{[m — (m — l)u m — a]a lm -(l-ijj m )a' im } cosrrnp rn + \ {m(l — w m )a 4m - [m(l— w m ) + l - a]a 2m } sin mip m ). 

where = w c sini and ip m = 0obs — 7r/2 — cp m (R). 



A2 Deprojection with incorrect viewing angles and incorrect centre 

Now we investigate how the expansion of the line-of-sight velocity vi os = Ci (R) cos iip + Si (R) sin iip changes if an incorrect 
inclination i, position angle F or centre (x, y) is chosen. We will derive a general equation for the line-of-sight velocity which 
is also valid if a combination of these incorrect choices occur, since to first order we can simply add the contributions of the 
individual incorrect parameters. 

Let us start with an incorrect inclination i = i+Si, where Si is the deviation from the correct value of the inclination i and is 
assumed to be small. Then q = cos(i+<5i) = q+8q. Using the equations for a projected circular orbit x" — Rcos tp, y" = qRsin ip 

r r/2 _|_ yiii jq2 we g n( j that the relation between cosip and cost/; is 



and cos?/> = x" / R — x" /yi 
cosip — (1 
Or, more general 
cos aip 
and 

fir, 

- cos 2ip)\ 



— ) cos ip H cos Sip. 

4q 4q 



aSq 



cos aip + (cos(a + 2)ij> — cos(a — 2)ip) + 0(Sq 2 ). 



fl = Aii + g(i 



(A14) 



(A15) 



(A16) 



Note that the amplitudes of the sine and c osine terms in equation (A12) are first order, except the cosip term whose amplitude 
is 1, and the corrections in equation (A15) are first order too. Therefore, only the deprojection of the cosip term will give rise 
to first order contributions to the line-of-sight velocity. 

Now define F = T + ST, where F is the correct position angle and ST the deviation from it. In the same way as for the 
inclination we find 



cos ip — cos ip - 
More general 



ST 
T 



(3q+-)sinip + (- 



q) sin Sip 



7 aST 
cos aip = cos atp H — - — 
Iq 

and 



R = R[l + i(i - q) sin 2^<Sr] 



i(a — 2)tp + sin(a + 2)^] - (l + q 2 )smaip L + 0(ST 2 ) 



(A17) 



(A18) 



(A19) 
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Again, only the deprojection of the cosip term will give rise to first order effects. 

Finally let us consider an incorrect choice of the kinematic centre. Suppose the centre of the coordinate system is shifted 
from (0, 0) to (x' , y'). Let us write this shift as (Sx, Sy) where the x-axis coincides with the major axis of the projected orbit. 
Then 



* SX A Si! A 

cost/; = cost/> + — — (1 — cos 2tp) — — — sin 2^. 
In general: 

cos axfj = cos aijj + [cos(a — l)ip — cos(a + 1)V>] + ^f^L [sin(a — 1)t/> — sin(a + 1)$1 + G(Sx 2 , Sy 2 ). 



2qR 



(A20) 



(A21) 



and 

R = R[l + ^-cosip + ^rSmip) (A22) 
R qR 

Here too only the deprojection of the cos rp term is important. 

Define dR = R — = — cos 2$) + j(~ — q)5T sin 2"ip + £| cos $ + ^4 sin^] and substitute the relations between tp 

and i/ 1 an d R and i? into the expansion of vi os . To first order a(R) = Ci(i? + dR) = d(R) + dR(dd(R) / dR) . Since da/dR is 
of the same order as d and dR is also small (order Sq, ST, 5x,8y), the only relevant term is dci/dR, the shape of the rotation 
curve (so this contribution is zero in the case of a flat rotation curve). Noting that Rdci(R) / 'dR — v*(R)a we can expand the 
line-of-sight velocity in the case of incorrect viewing angles and centre as: 



I Sx Sq « ST 



(Sq+-)-a(--q) 



sin ip — (1 — a) — — cos 2^> 
2R 



-(1 - a) ^ sin 2^ + (l-a)^cos3^ - (l-a)^(i - 9 )sin3^ [■ - V/^ '* <•«,- /V - *,('/>') mu I, : (A23) 



i>0 



If this is expanded into new harmonics, v\ OB = ^\ Cj(i?) cos ii/> + Sj(iZ) sini^>, then Sj — Si,Ci = a, except for i < 3. 
A3 Deprojection under the assumption of a circular velocity field 

Now we will calculate the errors in F, i and centre if the general velocity field is deprojected under the assumption of a circular 
velocity field ui os = w, costp. As a result, the x 2 deviation from the best fitting circular field is 



X 2 (Sq, ST, Sx, Sy) ~ ^(<= 2 + sf) - c{, 

i>0 



(A24) 



where §i and c; are expressed in terms of the Sq, ST, Sx and Sy in equation (A.23). The best fitting values of ST and Sq are 
given by 



Sq=- 



4qc 3 



(ci +c 3 )(l 



<5F = 4g 



[(3g 2 + l)-a(l-g 2 )] si + (1 - g 2 )(l - a) 



S'3 



{[3g 2 + 1 - (1 - g 2 )a] 2 + (1 - «) 2 (1 - g 2 ) 2 ]} (ci + c 3 ) ' 



Equation (A.24) is minimized with respect to Sx and Sy when 

6x = 2Rc 2 /{l-a), 6y = 2qRs 2 /(l-a). 

The resulting expressions for d, Si under a circular fit are: 



1 + a 

Co = Co + C2, 



c 2 =0, 



1-a 



<5<Z, 



§2 =0, 



ci = ci [1-^2.(1 -a)] = ci +c 3 , si =si - ci [(3g + - (1/q - q)a] 5T/4, 



£3 = c 3 + (1 - a)^ci = 0, 
4g 



s 3 = S3-ci(l/g-5)(l-a)jr/4. 



(A25) 



(A26) 



(A27) 



j43. i Effect of m — 1 distortion 



For an m = 1 distortion equation (A13) gives the correct coefficients if the correct centre is chosen: 



co = v* \ [(1 — 2oji + a)an + (1 — a>i)a 3 i] cos 931 + |[(1 — 0^1)041 + (a — wi)a2i] sin</>i, 



C2 
S2 



y*3 [(1 - a)an - (1 — ^i)a 3 i] cosi^i + |[(1 — ^i)a 4 i - (2 — U\ — a)a 2 i] sini^i 



v*| [(1 — a)on - (1 - wi)o3i] sini^i + |[(1 — wi)a4i — (2 — uii 



a)a2i\ cos (pi. 



(A28) 
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If the velocity field itself is used to derive the centre, these equations reduce according to equation (JA27|) to: 
s ° = v * 2(e»-i) {[(1 — «)^ion - (1 - wj)o3i] cosy>i — 2[(1 - wi)o4i + («wi - l)oai] sin^i} , 

c 2 = 0, (A29) 

S 2 = 0. 

These terms are not affected by an error in q and V. 
A3. 2 Effect of m = 2 distortion 



With m = 2, equation (A13) gives the coefficients if a correct inclination and position angle are chosen: 
ci = v* {l — | [(2 - 3o2 + a)ai2 + (1 - 02)032] cos 2<^2 + § [2(1 - 02)042 + (1 — 202 + 0)022] sin 2<^2 } , 
si = «* { — j [(2 — 3o 2 + a)ai2 + (1 - 02)032] sin2^ 2 + \ [2(1 — 102)042 + (1 - 2o 2 + 0)022] cos2<^ 2 } , 

C3 = n» { — j [(2 — CJ2 — a)«i2 — (1 — ^2)032] cos 2(^2 + ^ [2(1 — 012)042 — (3 — 2^2 — 0)022] sin2ip2} , (A30) 
s 3 = v, { — j [(-2 + o 2 + a)ai2 + (1 - 02)032] sin 2952 + \ [2(1 - 02)042 - (3 - 2o 2 — 0)022] cos2<y3 2 } . 



sin 2<p 2 



When the velocity field itself is used to derive the position angle and inclination, we find using equation flA27p 

£1 = w» {l + I [(a — 02)012 + (1 - 02)032] cos 2^2 + [2(1 - 02)042 + (a - 1)022] sin2<p 2 } , 

si = w* { — \ [(2 — 3o 2 + o)ai2 + (1 - 02)032] sin 2(^2 + 5 [2(1 — 02)042 + (1 — 2o 2 + 0)022] cos 2^2 — 

[3o 2 +l+a(g 2 -l)]({2(l-Q)(o 2 -o) + 2[2-a 2 -4o 2 +a(2-H.;2)]g 2 } oi2-[l + g 2 -a(l-g 2 )](l^ 2 )o 3 2) 
4[1 + 2g 2 + 5g 4 + a 2 (l - g 2 ) 2 + 2a(g 4 - 1)] 

[3g 2 +l+a(g 2 -l)] {4[l+g 2 -a(l^g 2 )](l^ 2 )a 42 - 2[l-(3-4o 2 )g 2 -a(l+g 2 )]a 2 2} 1 
4[1 + 2g 2 + 5g 4 + a 2 (l - g 2 ) 2 + 2a(g 4 - 1)] C ° S 2v " 2 J ' 

£ 3 = 0, (A31) 
h = w* { — \ [(02 + a - 2)012 + (1 - 02)032] sin 2(^2 + 5 [2(1 — 02)042 — (3 — 2o 2 - 0)022] cos2</3 2 + 

(l-a)(l-g 2 )({2(l-a)(o 2 -Q) + 2[2-a 2 -4o2+a(2-k.;2)]g 2 } ai2-[l+g 2 -a(l-g 2 )](l-o;2)a32) 
+ 4[1 + 2g 2 + 5g 4 + a 2 (l - g 2 ) 2 + 2a(g 4 - 1)] Sm 2 ^ 2 + 

(l-a)(l-g 2 ) {4[l^ 2 -a(l-o 2 )](l-^2)a 4 2-2[l-(3-4o 2 )g 2 -a(l+g 2 )]o22} 1 
+ 4[1 + 2g 2 + 5g 4 + a 2 (l - g 2 ) 2 + 2a(g 4 - 1)] C ° S ^ J ' 

These equations are independent of a small error in the position of the centre. For a triaxial non-rotating halo model, in which 
the potential is scale free, we have 02 = 0, a = (flat rotation curve), the phase of the perturbation is constant with radius 
<fi2(R) ~ const (therefore 042 = 022 = 0), 032 = 2ai2 = 2e pot (FvGdZ) and these equations simplify to the equations ([]). 

A3. 3 Solid body rotation 



Let us finally consider the effect of solid body rotation, a = 1. Equations (A.25) and ( |A2q ) show that in this case it is 
impossible to measure the inclination. Indeed, in solid body rotation, all iso-velocity contours are parallel and each inclination 
fits equally well. The position angle on the other hand can be determined, ST = Centre and systemic velocity cannot be 
determined uniquely. Therefore, assumptions for all these parameters have to be made in order to get a rotation velocity. It 
is not possible to interpret harmonic terms under these conditions in the way we did in this Appendix. 
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